Projection Operator Approach to 
Spin Diffusion in the Anisotropic Heisenberg Chain at High Temperatures 
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We investigate spin transport in the anisotropic Heisenberg chain in the limit of high temperatures 
(/3 — !■ 0). We particularly focus on diffusion and the quantitative evaluation of diffusion constants 
from current autocorrelations as a function of the anisotropy parameter A and the spin quantum 
number s. Our approach is essentially based on an application of the time-convolutionless (TCL) 
projection operator technique. Within this perturbative approach the projection onto the current 
yields the decay of autocorrelations to lowest order of A. The resulting diffusion constants scale as 
A~^ in the Markovian regime /S. ^ 1 {s — 1/2) and as A~^ in the highly non-Markovian regime 
above A ~ 1 (arbitrary s). In the latter regime the dependence on s appears approximately as an 
overall scaling factor s{s + 1) only. These results are in remarkably good agreement with diffusion 
constants for A > 1 which are obtained directly from the exact diagonalization of autocorrelations 
or have been obtained from non-equilibrium bath scenarios. 

PACS numbers: 05.60.Gg, 05.30.-d, 05.70.Ln 



Transport in one-dimensional quantum systems has 
been a topic of theoretical investigations since at least 
the 1960s P, Q- Even nowadays there is an ongoing and 
still increasing interest in understanding the transport 
phenomena in such systems, including their dependence 
on length scales and temperature [3l-[l5| . Much work has 
been devoted to a rather qualitative classification of the 
emerging transport types into either non-normal ballistic 
or normal diffusive behavior. In this context the crucial 
mechanisms for the emergence of pure diffusion have been 
intensively studied and in particular non-integrability or 
quantum chaos are frequently discussed w.r.t. their role 
as an at least necessary prerequisite [l2-14|. Spin chains 
are a central issue of research with a considerable focus 
on spin transport in the anisotropic Heisenberg chain, in 
particular on the case s = 1/2 |3l4l3l|. In that case there 
still is an unsettled debate about the concrete range of 
anisotropics where the dynamics at finite temperatures 
is indeed diffusive 3-5]. 

However, even if the diffusive range of anisotropics was 
definitely known, the question about the quantitative 
value of the diffusion coefficient arises naturally. In fact, 
this question is as challenging as the proof of diffusion as 
such. Because the majority of all available approaches to 
transport has to deal more or less with finite quantum 
systems, signatures of diffusion in the thermodynamic 
limit may not be observed, e.g., due to large mean free 
paths, even in the limit of high temperatures. For these 
temperatures diffusion coefficients can be found in the 
literature for s = 1 and the isotropic point [H 0, [3l and 
for s — 1/2 and a narrow window of large anisotropics 
above the isotropic point [7-ll| solely. On that account 
there is an urgent need for a much more comprehensive 
picture of diffusion constants. 



In this work we intend to make an essential step towards 
such a picture at high temperatures. To this end we focus 
on the quantitative evaluation of diffusion constants from 
current autocorrelations, either directly or perturbatively 
by the use of projection operator techniques, see Refs. 
Ig- 1^ Within this perturbative approach the projection 
onto the current yields the decay of autocorrelations to 
lowest order of the anisotropy parameter. The resulting 
diffusion constant is a smooth function of the anisotropy 
and agrees well with the above ones from the literature 
for large anisotropics, i.e., at the upper boundary of the 
expected range of validity. This agreement also supports 
the validity of the perturbative approach in the limit of 
small anisotropics, if transport is indeed diffusive. 

The anisotropic Heisenberg chain (XXZ model) may be 
described by a Hamiltonian of the form H = Hq -|- A y, 
where Hq and V are given by 
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N 



M+1 



(1) 



with periodic boundary conditions. Here, J denotes the 
nearest neighbor coupling strength, A the anisotropy, N 
the number of sites, and the operators {i = x, y, z) 
represent the standard spin matrices (w.r.t. site fj.) for the 
spin quantum number s. For the clarity of later notations 
we define the abbreviation s = •y/s(s + 1). 

The Einstein relation (Tdc = X ^ connects the spin (dc) 
conductivity (Tdc and the spin diffusion constant T) via 
the static susceptibility x- Therefore the spin diffusion 
constant according to linear response theory (LR) 
reads in the limit of high temperatures (/3 0) 



V = lim V{t) , 



X dimH N 



dt'C{t'), (2) 



2 



where x = 1/3 /3s^, dim^ = (2s + 1)^, and the spin 
current autocorrelation function is given by 



C{t) = Tr{ J(t) J(0)} , J 



N 



with the initial value Tr{j2} = 2/9 J"^ diuYHN . 



(3) 



For 



s = 1/2 the infinite time integral in Eq. ^ diverges for 
A = due to [Hq, J] = [2l[. It also diverges, whenever 
there is only partial conservation, i.e., a non-zero Drude 
weight. Finite Drude weights indicate ballistic transport 
at the infinite time scale but are commonly expected to 
vanish in the thermodynamic limit for the non-integrable 
cases s > 1/2 [l^. For the integrable case s = 1/2 Drude 
weights are widely expected to vanish for A > 1 d, 0| , 
while they may already be zero for A = 1 0, Q] or even 
for ah < A < 1 [i]. 

Apart from the common LR picture, the time-dependent 
diffusion coefficient V{t) according to Eq. © can be also 
connected to the time evolution of the spatial variance of 
an initially inhomogeneous non-equilibrium density, see 
Refs. [13, El Within this picture the density dynamics 
is diffusive at a certain time, respectively length scale, if 
is constant at that scale. Such a scale may exist, 
even if 2?(i) eventually diverges in the infinite time limit, 
i.e., despite finite Drude weights. We thus do not focus 
only on the infinite time limit, but investigate the full 
time dependence of X'(t), too. 
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FIG. 1: The time dependence of the spin diffusion constant 
"Dit) in the limit of high temperatures (/3 — >■ 0), as obtained 
directly from exact diagonalization. The data is displayed for 
the anisotropy parameter A = 1.5 and for the spin quantum 
numbers s = 1/2 (iV = 20), s = 1 (iV = 12), s = 3/2 (N = 9), 
and s = 2 (A'^ = 8). The data in the inset is shown for the 
spin quantum number s = 1/2 and for various chain lengths 
iV = 8, 10, ...,20. 

A possible strategy for the analysis of the dependence 
of 'D(t) on time is the direct application of numerically 
exact diagonalization (ED). Of course, ED is restricted 



to a rather limited range of A'', even if all symmetries are 
taken into account, e.g., translation, rotation about the 
z-axis [i^l . But finite N effects are often less pronounced 
at finite time scales. In Fig. [T]we show 'D{t) at such time 
scales for A = 1.5 and the maximum N which is available 
to us for a certain s: 'D{t) firstly increases at short times 
t < 2/(Js) and then becomes approximately constant 
at intermediate times, i.e., a "plateau" with a height of 
about V w 0.7 J s is formed here. 2?(t) finally increases 
again due to non-zero Drude weights for finite TV. But in 
the thermodynamic limit Drude weights are expected to 
vanish for A = 1.5 and each s. Furthermore, the height 
of the plateau does not change with N, while its width 
gradually increases, see Fig. [T] (inset). The height of this 
plateau consequently is a reasonable suggestion for the 



concrete value of the diffusion coefficient 
for s = 1/2 the value 2? « 0.6 J is in remarkably 



Mim. In fact, 

„ood 

agreement with non-equilibrium bath scenarios [8|, |9| . 
Analogously, the concept of plateaus allows to suggest 
further values of the diffusion coefficient for anisotropy 
parameters above A 1.4, see Fig. [3] (squares). Here, the 
dependence of 2? on s appears approximately as a scaling 
factor s. For anisotropy parameters below A ~ 1.4 the 
suggestion of diffusion coefficients is not possible, because 
Drude weights become dominant and plateaus disappear 
for finite N. Plateaus appear only, if contributions from 
the Drude weight and the regular part are well-balanced, 
at least to some degree. However, plateaus are also not 
visible, if both contributions are treated separately from 
each other ill]. 

Another strategy for the analysis of the time evolution 
of C{t), respectively T>{t) is provided by an application 
of the time-convolutionless (TCL) projection operator 
technique [l6l . [l7| . This technique, and the well-known 
Nakajima-Zwanzig (NZ) method 18, 13], are commonly 
applied in order to describe the reduced dynamics for 
a set of relevant variables. Their application essentially 
requires a Hamiltonian of the form H — Ho + AV and 
the commutation of the relevant variables with Hq (or a 
comparatively slow dynamics w.r.t. Hq). Although both 
methods are well-established approaches in the context of 
open quantum systems, the following approach to current 
autocorrelation functions in closed quantum systems is 
an uncommon concept. However, in this context NZ is 
very similar to the more common Mori-Zwanzig memory 
matrix formalism, see Ref.0- But a "Mori-TCL" variant 
does not exist, of course. 

The above-mentioned set of relevant variables is specified 
by the definition of a suitable projection (super)operator 
V which projects out the relevant part of an operator 
p{t). (This operator does not need to be a density matrix 
in the strict sense.) To this end we define 

A^TrlPyU, Ppit)^Tr{Am}A, (4) 

where the normalized current operator A is introduced 
in order to satisfy the property = "P of a projection 



3 



(super)operator. 

For the initial condition p(0) = A the TCL technique 
routinely yields a homogenous differential equation for 
the actual expectation value a{t) = Tr{A{t) Ai{t)}. Here, 
the index I of Ai(t) denotes the interaction picture, i.e., 
the Heisenberg picture w.r.t. Hq. Since C{t) oc a(t) for 
Ai{t) — A (or Ai{t) « v4 at a pertinent time scale [2l|). 
the TCL equation is identical for C{t) in that case. It 
then reads 



dt 



i=l 



and avoids the often troublesome time-convolution which 
appears in the NZ variant [l6, 17|. The time-dependent 
rate R{t) is given in terms of a systematic perturbation 
expansion in powers of A. (All odd orders of A vanish for 
this and many other quantum systems.) The truncation 
of R{t) to lowest order of A reads 



R{t)^A^ f dt'fd') 
Jo 

with the two-point correlation function 

fit) = Tr{ J2}-1 Tr{ *[J, VUt) i[J , l/]i(0) } , 



(6) 



(7) 



where the above index I indicates again the Heisenberg 
picture w.r.t. Hq. In general this lowest order truncation 
is expected to be justified for small A (and short t). Even 
though the incorporation of higher order corrections is in 
principle possible, their concrete evaluation is an almost 
impossible task, both from an analytical and numerical 
point of view. 

The concrete evaluation of the lowest order truncation for 
R{t) in Eq. ([6|) is feasible. Moreover, this evaluation may 
be done analytically, if Hq can be brought in diagonal 
form, e.g., via the Jordan- Wigner transformation onto 
non-interacting spinless fermions for s = 1/2. However, 
since the expression for R{t) still takes on a non-trivial 
form after such a transformation, it may possibly have 
to be evaluated numerically at the end. We thus directly 
evaluate R{t) in Fig. [2] numerically by the use of ED, 
similarly as done in Refs. Q 

There is an apparent similarity between Figs. [Hand [21 
even though completely different quantities are shown. In 
particular a constant rate R « 0.67 A^ J s may be read 
off from the height of the plateau at intermediate times, 
at least for s — 1/2. Of course, for s > 1/2 the plateau is 
not developed as clearly for finite iV. But the respective 
i?(t)-curves have converged completely for short times 
t < 2/{Js) solely. We hence assume a constant rate R 
in that case also, e.g., as given by the height of the local 
maximum at such times, see Fig. [2| (arrow). Obviously, 
this assumption implies a positive R, i.e., we expect that 
TCL does simply not predict Drude weights to lowest 
order, see below. However, in principle R{t) may become 
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FIG. 2: The time dependence of the decay rate R{t), as given 
by Eq. ([6]). The data is shown for the spin quantum numbers 
s = 1/2 (iV = 20), s ^ 1 [N = 12), s = 3/2 (N = 9), and 
s = 2 {N = 8). The data in the inset is displayed for the 
spin quantum number s = 1/2 and for various chain lengths 
iV = 8, 10, . . . , 20. The extracted plateau height for s > 1/2 
is indicated (arrow). 



zero in the long-time limit. 

The above assumption is relevant only in the Markovian 
limit of long relaxation times, i.e., small A. In that limit 
the rate R{t) can be replaced by the constant value R and 
lowest order TCL consequently predicts the exponential 
decay C{t) = C(0) exp(-i?t), cf. Eq. ©. The diffusion 
constant in Eq. ^ eventually scales as P oc R~^ oc A^^, 
e.g., as generally expected from a Boltzmann equation 
approach to the same question. The Markovian limit of 
long relaxation times concretely turns out to be realized 
for A < 1, see Fig. [3] 

However, in that limit the resulting diffusion constants 
for s > 1/2 have to be understood in terms of an upper 
boundary. Since Ai{t) « A does not hold true for long 
relaxation times, a back-transform from the interaction 
picture becomes necessary in order to obtain C(t) from 
a(t). Because we do not perform this back-transform due 
to its unavailability, we only consider the contribution 
from V and neglect any contribution from Hq, i.e., C{t) 
probably decays too slowly. Obviously, C{t) decays too 
slowly, if exclusively Hq contributes, i.e., for A = 0. 
For small A there is no substantial difference between 
the predictions of lowest order TCL and NZ. But the 
situation changes for large A, when relaxation times are 
short and non-Mar kovian effects become relevant, i.e., 
the initial increase of the rate R{t) at short times, see 
Fig. [2] At such times the R{t)-cmve has well converged 
already for each s and additional assumptions like the one 
above are not required here. Since the initial increase 
appears to be approximately linear, lowest order TCL 
predicts a Gaussian decay in the highly non-Markovian 
limit. As a consequence the diffusion constant scales as 
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FIG. 3: The spin diffusion constant as a function of the 
anisotropy parameter A in the limit of high temperatures 
(/? — >■ 0), as predicted by TCL to lowest order of A, for the 
spin quantum numbers s = 1/2, . . . , 2 (curves). Further data 
from exact diagonalization (squares, s > 1/2) and from the 
literature are indicated: Ref. |3i,rhombus, s = 1/2), Refs.|l,i 
(triangles, s — 1/2), and RefsJU. [^. [TH (circle, s = 1). Curves 
for s > 1/2 and A 1 have to be understood in terms of 
an upper boundary. {D is located somewhere below, e.g., as 
sketched by the filled area.) 



D oc A^^ now. The highly non-Markovian limit of very 
short relaxation times turns out to be concretely realized 
for anisotropics above A ~ 1, see Fig. [31 

For s > 1/2 the diffusion constants according to lowest 
order TCL in Fig. |3] are in excellent agreement with the 
results from ED for A > 1 in Fig. [T] Furthermore, these 
diffusion constants are in good agreement with results 
for A = 1 in the literature (standard frequency moments 
analysis, microcanonical Lanczos method P, 0, For 
3 = 1/2 there is a deviation on the order of 20% between 
diffusion coefficients according to lowest order TCL and 
results for A ~ 1.5 from ED in Fig. [T] or the literature 
(non-equilibrium bath scenarios [1, @| ) ■ Such a deviation 
probably arises from the neglection of all higher order 
corrections which are more important for s = 1/2. 
Of course, the above agreement also supports the validity 
of lowest order TCL for s = 1/2 and A <C 1, e.g., the 
usual range of application for a perturbation theory. But 
one has to keep in mind that this perturbation theory is 
generally restricted to small A and short t as well. The 
restriction becomes manifest, since we expect that TCL 
does simply not predict Drude weights to lowest order, 
i.e., the non-decaying character of C{t) essentially is a 
higher order effect, as known for large A. Consequently, 
lowest order TCL can be valid at all t only, if Drude 
weights eventually vanish in the thermodynamic limit. 
On that account we may finally conclude as follows: 
Whenever there is spin diffusion in the XXZ model at 
high temperatures, either at a finite or the infinite time 
scale, we expect that the respective diffusion coefficient 



is reasonably suggested by lowest order TCL, e.g., as 
displayed for s = 1/2 in Fig. |3l 
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